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SUMMARY 


Formulae  and  prooesses  are  presented  for  determining  the  beat  shape  for 
the  external  surface  of  an  annular  duct  in  order  to  produce  a minimum  amount 
of  wave  drag  in  supersonic  flow  under  certain  conditions  which  are  invoked  in 
order  to  ensure  that  the  missile  geometry  will  obey  practical  design  require- 
ments. The  question  of  what  the  best  shape  is,  at  least  to  the  degree  of  ap- 
proximation inherent  in  the  use  of  the  linear  theory,  is  solved  under  the 
restriction,  first,  that  the  area  between  the  (given)  inner  contour  of  the 
body  of  revolution  and  the  sought  outer  contour  is  a constant,  and,  second, 
under  the  stipulation  that  the  volume  conprised  between  the  surfaces  whioh 
cure  swept  out  when  such  inner  and  outer  contours  are  rotated  about  the  duct 
axis  is  to  be  a constant.  Two  further  distinctions  are  also  made  in  the 
treatment  acoorded  this  problem  on  the  basis  of  the  type  of  inner  (known)  duct 
shape  given;  i.e«,  in  one  case  it  is  assumed  that  the  annular  duct  will  dif- 
fer but  slightly  from  a cylinder  while  in  the  other  case  it  is  assumed  that 
the  basic  shape  upon  which  the  desired  external  duct  oontour  is  to  be  built 
up  is  fundamentally  a frustum  of  a cone.  The  flow  through  the  inside  of 
the  duot  is  not  considered  at  all,  and  the  internal  oblique  shock  is  thus 
assumed  to  be  attached  to  the  entrance  lip;  the  external  contour  is  also 
assumed  to  *olose*;  i.e.,  the  trailing  edge  at  the  duct  exit  does  not  have 
a blunt  face.  This  problem  had  previously  been  solved  when  no  perimetral 
condition  was  invoked,  but  under  the  more  stringent  conditions  now  being  im- 
posed it  is  believed  that  more  directly  useful  results  will  be  forthcoming 
when  the  numerical  applications  are  made  in  another  report  to  follow.  Two 
distinct  modes  of  attack  are  offered  to  attain  the  solution  in  all  oases; 


1.  INTRODUCTION 


The  objeot  of  this  study  is  to  determine  what  is  the  best  shape  for  an  annular 
duot  when  it  is  to  operate  in  a supersonic  stream  under  the  usual  stipulation  that 
linear  theory  is  sufficiently  aoourate  for  the  purposes  in  mind.  In  other  words,  the 
problem  resolves  itself  into  finding  out  what  form  the  meridional  line  which  delimits 
the  external  oontour  of  the  annular  duct  in  question  must  have  in  order  that  the 
external  drag  due  to  this  duct  shape  should  be  rendered  a minimum,  under  either  of 
the  following  two  perimetral  conditions} 

(a)  the  area  enclosed  between  the  external  meridional  line  defining  the  outside  oontour 
of  the  duct  and  the  internal  meridional  line,  whioh  is  assumed  to  be  known  beforehand, 
is  required  to  have  a set  value. 

(b)  the  volume,  enclosed  between  the  surfaces  of  revolution  swept  out  when  the  above 
oited  meridional  lines  are  rotated  about  the  duet's  axis,  is  required  to  have  a oertain 
fixed  value. 

The  problem  will  be  attacked  in  two  stages.  First,  the  case  where  the  radius  of 
the  circular  cross-section  (or  parallel)  representing  the  entrance  of  the  duct  is  equal 
to  the  radius  of  the  parallel  representing  the  exit  of  the  duot  will  be  treated  (see 
sketoh  *a"  in  Fig.  1),  and  then,  after  that,  the  more  general  case  will  be  examined, 
wherein  the  entranoe  and  exit  cross-sections  may  have  unequal  radii  (see  sketch  *b*  in 
Fig.  1). 


on*  la  baaed  on  the  uae  of  Ll(hthlll'a  V-funotlon  and  givoa  the  boat  oontour 
ahapa  aad  drag  dlreotly,  while  tha  alternative  method,  a tanning  from  pra- 
vlouily  publlabad  lnveatlgatlone  of  tha  author,  raqulraa  that  a aultabla  dia- 
trlbutlon  of  auparaonlo  aouroaa  ba  dataminad  firat  of  all,  aad  than  onoa  tha 
daaorlptlon  of  how  thaaa  aingularltiea  vary  along  tha  x-axla  la  obtalnad,  by 
a relatively  alnpla  prooaaa,  tha  baat  oontour  and  tha  ralatad  drag  are  derived 
aa  auxiliary  infornatlon.  Tha  oonneotion  between  thoae  raaulta  and  thoaa  whioh 
would  ba  obtained  by  uaa  of  two-dlmenaional  (Aokeret)  theory  la  alao  pointed 
out  aa  obiter  dlota. 
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duot  when  it  is  to  operate  in  a supersonic  stream  under  the  usual  stipulation  that 
linear  theory  is  sufficiently  accurate  for  the  purposes  in  mind.  In  other  words,  the 
problem  resolves  itself  into  finding  out  what  form  the  meridional  line  which  delimits 
the  external  contour  of  the  annular  duct  in  question  must  have  in  order  that  the 
external  drag  due  to  this  duot  shape  should  be  rendered  a minimum,  under  either  of 
the  following  two  perimetral  conditions  s 
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is  required  to  have  a set  value. 

(b)  the  volume,  enclosed  between  the  surfaces  of  revolution  swept  out  when  the  above 
oited  meridional  lines  are  rotated  about  the  duet's  axis,  is  required  to  have  a certain 
fixed  value. 

The  problem  will  be  attacked  in  two  stages.  First,  the  case  where  the  radius  of 
the  ciroular  oross-seotion  (or  parallel)  representing  the  entrance  of  the  duct  is  equal 
to  the  radius  of  the  parallel  representing  the  exit  of  the  duot  will  be  treated  (see 
sketoh  *a*  in  Fig.  1),  and  then,  after  that,  the  more  general  case  will  be  examined, 
wherein  the  entranoe  and  exit  cross-sections  may  have  unequal  radii  (see  sketch  *b*  in 
Fig.  1). 
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Coordinate  Syetem  and  Geometric  Relatione  Applying  to  the  Caee  Where  the  8ought 

Annular  Duct  ie  Cylinder-Like  in  8hape 

Pig.  1(a) 

kR 
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Coordinate  8yetem  and  Geometric  Relatione  Applying  to  the  Caee  Where  the  Sought 
Annular  Duot  ie  Cone-Like  in  3hape 

Pig.  1(b) 
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Two  different  approaohes  will  be  followed  in  obtaining  the  solution  to  the 
problems  thus  set.  One  will  be  based  on  the  application  and  extension  of  the 
work  done  by  Lighthill  in  Referenoe  1 (this  way  of  handling  the  problem  will  con- 
stitute Part  I of  this  study) , while  the  other  approaoh  will  be  founded  on  the 
applloation  and  extension  of  the  investigations  oarried  out  by  the  author  pre- 
viously in  Referenoe  2 (this  treatment  will  oonstitute  what  is  oalled  Part  II  of 
the  present  investigation).  By  use  of  the  first  of  these  two  avenues  of  approaoh 
one  finds  what  the  external  contour  shape  has  to  be,  to  give  minimum  drag,  di- 
reotly,  through  solution  of  a non-homogene ous  integral  equation  of  the  seoond 
kind.  In  the  oase  of  the  seoond  proposed  way  of  looking  at  the  problem,  however, 
one  has  to  determine  first  of  all  what  the  distribution  of  singularities  (super- 
sonio  souroes)  along  the  axis  of  the  duot  has  to  be,  and  it  is  found  that  the 
description  of  this  source  distribution  is  given  by  solution  of  an  integral  equa- 
tion of  the  first  kind. 
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2.  LIST  OF  PRINCIPAL  SYMBOLS 


distance  along  tho  ixli  of  the  duot,  measured  positive  In  the  dlreotion 
of  tho  free  stream 

froo  itream  velooity  vootor 

Porpondloular  dlotanoo  away  from  tho  duot  axil 

radius  of  tho  inlot  aootion  of  tho  duot 


radius  of  tho  outlet  oootien  of  the  duot 
*4  * Rg 
2 

radiua  of  the  Inner  oiroular  aootion,  produood  when  tho  annular  duot 
ia  out  by  tho  plane  x ■ oonatont 

radiua  of  tho  outer  oiroular  aootion,  produood  when  the  annular  duot 
ia  out  by  tho  aamo  plane  x • oonatant 

Rg  m R]_ 

E 


Uaoh  number  of  tho  froo  stream  flow  with  volooity  U_ 

X m X 

*1'  *2  “ 

BRj  B Rg 

*•  * ri  • , r2  - *• 

; 1 TTIT  2 


length  of  tho  duot,  measured  along  its  axis 

Jl  .<i  ■ _£_.•*!  -JL 

BR0  BRi  BRg 

dri  • dRg 
da  dx 

looal  strengths  of  tho  aupersonlo  aouroe  distribution  whloh  la  plaood 
along  tho  axis  of  tho  duot,  and  whloh  is  responsible  for  tho  oroation 
of  tho  porturbod  flow  whloh  flews  past  and  in  oonformity  with  tho 
external  oontour  of  tho  annular  duot 


LIST  or  PKIHOIPAL  SYMBOLS  - (Oont.) 


P 

Pm 


r 


■ local  proaiuro 

* proa aura  oorraapondlng  to  a point  looatad  in  tha  vmiform  undiaturbad 
fraa-atraaa  flow 

* looal  danalty 

* danalty  oorraapondlng  to  a point  looatad  In  tha  unlfom  undiaturbad 
fraa-atraam  flow 

■ adlabatlo  axponantj  aquala  1.406  for  air 

■ drag  ooafflolant  » drag  axartad  on  tha  axtarlor  aurfaoa  of  tha  duot 

i/>~  ui  R02 


3.  Cm  of  an  Annular  Duct  with  Cylindrical  Type  of  Central  Channel 
(Rl  * R2  * Re).  Solution  by  Bxtonfllon  of  Llghthlll^  Method 


The  drag  oooffiolent  that  reeulte  from  determination  of  the  sum  of  all  the 
resultant  foroes  oreatod  by  aotion  of  the  pressure  Increments,  p - p„f  that  are 
exerted  against  the  external  surfaoe  of  thi3  sort  of  duot  has  been  obtained  by 
Llghthill  In  Referenoe  1 and  the  pertinent  formula  has  the  form 


7?2  -2 


l/V  C U - ^/)  ?j>  ( *4 ) cf-d  d 


where  the  funotion  W(y) , as  a variable  that  is  dependent  on  the 
y-values,  has  been  evaluated  by  the  Admiralty  Computing  Service, 
lated  results  have  been  given  in  the  referenoe  by  Lighthill  just 
they  are  reproduced  here  in  the  appended  Table  I, 


? (1) 

assigned 

and  tho  tabu- 
mentioned  i 


TABLE  I 


Lighthill 's  W(y)  Funotion 


y 

mmm 

y 

w(y) 

0.0 

0.50000 

0.06929 

0.2 

0.43190 

0.06384 

0.4 

0.37552 

0.05896 

0.6 

0.32889 

BSS 

0.05068 

0.8 

0.28867 

4.8 

0.04370 

1.0 

0.25497 

5.2 

0.03803 

1.2 

0.22617 

5.6 

0.03331 

1.4 

0.20143 

6.0 

0.02936 

1.6 

0.18009 

6.4 

0.02602 

1.6 

0.16156 

6.8 

0.02320 

2.0 

0.14542 

7.2 

0.02078 

2.2 

0.13131 

7.8 

0.01869 

2.4 

0.11889 

8.0 

0.01688 

2.6 

0.10798 

8.4 

0.01534 

2.6 

0.09833 

8.8 

0.01397 

3.0 

0.08976 

9.2 

0.01277 

3.2 

0.08216 

9.6 

0.01170 

3.4 

0.07536 

10.0 

0.01077 

Let  attention  now  be  turned  to  the  perimetral  oondition  that  may  be  stated 
ae  follows:  in  any  arbitrary  general  one  of  the  meridional  half -planes  passing 


through  the  axis  of  the  body  of  revolution  under  discussion,  the  area  that  is 
to  be  enolosed  between  the  meridional  line  delimiting  the  outer  boundary  of  the 
annular  body  and  the  meridional  line  marking  out  the  inner  boundary  of  the  duot 
is  to  turn  out  to  have  a given  fixed  value.  This  oondition  is  expressible  by 
the  simple  relation 

~ B R.  0 + L 

( R.g  - R.  £ ) d x ~ constant 


On  the  other  hand  if  it  is  taken  for  granted  that  ■ R^(x)  is  known, 
then  this  constraining  oondition  will  have  the  form 


If  C -4)  Q ^ = C 


(2) 


O v0  VQ 

while  it  is  required,  in  addition,  because  of  the  assumption  that 


Rl  ■ R2  • R0,  that 


^ c 2 > - 0 


(3) 


must  hold  true. 

Looking  back  at  Eq.  (1)  now,  it  is  seen  that  when  an  arbitrary  variation 
A if  is  given  to  the  variable  the  first  variation  in  Cq  (whioh  will  be  de- 
noted here  by  the  symbol  A Cp)  will  take  on  the  form 

acd  = 

provided  one  uses  the  shorthand  notation  that 


f*  -4F  <■$)]  Ar,  <V 


(4) 


F*V  = 


/ 

(u  - ^ / ; 7 u ) 


(5) 
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Moreover,  sinoe  conditions  (2)  and  (3)  must  always  be  satisfied. 


the  var- 


iation in  , denoted  by  A rf  , must  obey  the  restrictions 

* 


and 


f A V C?)  m° 

l'*L 


(6) 


A (*)  c/4. 

If  the  order  of  the  integrations  indioated  in  the  last  of  the  restraining 
conditions  given  as  Bqs.  (6)  is  now  inverted,  and  if  due  advantage  is  taken  of 
the  relationship  set  down  as  the  first  of  the  Eqs*  (6),  then  the  perimetral  con- 
dition beoomes  just 

/ & T(  (^  ) 4*0  ~ O 

/ (?) 

Jo 

If  the  slope  function  ^ (z)  really  constitutes  the  description  of 
the  meridional  line  contour  which  produces  a minimum  drag  coefficient,  Cp,  then 
the  first  variation  A Cp  has  to  turn  out  to  be  zero  whatever  the  variation 

of  , A 7f  , is,  so  long  as  it  satisfies  the  conditions  stated  mathematically 
by  means  of  Eqs.  (6)  and  (7).  Consequently,  if  A ^ and  A 2 are  taken  to  be  two 
quite  arbitrary  oonBtants,  a priori,  it  follows  that  the  relation  given  below 
must  be  satisfied: 

8 7}  — 4 F (y)  + A / + ^ « O 

(8) 

If,  on  the  other  hand,  the  perimetral  condition  imposed  on  the  shape  of  the 

external  duot  oontour  being  sought  is  the  one  described  in  the  (b)  part  of  the 

first  paragraph  in  the  Introduction,  this  situation  is  translated  into  analytie 

language  by  writing  that 

8R0t  L CBRo  + L 

(R-l  - R})  c ix  -Z 

6R.  o 

must  be  invoked. 


(Re-  R't)  cfx  = constant 
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Now,  since  in  the  case  under  consideration  one  has  approximately  that 


P e +_  ^ i 

ZP-o 


1 


(9) 


it  is  seen  that  once  again  the  relation 

Jl 

(st  d % * constqnt  (2 * ) 

holds,  which  is  identical  with  condition  (2). 

It  is  evident,  therefore,  that  to  the  degree  of  approximation  premised  by 
writing  the  simplification  given  as  Eq.  (9),  the  perimetral  conditions  stated 
in  the  Introduction  aB  stipulations  (a)  and  (b)  are  both  represented  mathemati- 
cally by  the  same  equation,  to  all  intents  and  purposes  of  the  study  now  to  be 
undertaken.  Of  course,  it  is  true  that  there  does  exist  a difference  between 

the  aotual  constants  which  appear  on  the  right  hand  sides  of  Eqs,  (2)  and  (2*), 

\ 

but  there  is,  none-the-less,  no  difference  in  the  general  form  of  these  two 
equations  of  restraint. 

It  is  worth  while  to  reiterate  that  Eq.  (8)  delineates  precisely  the  meri- 
dional line  corresponding  to  a body  which  has  least  drag,  is  a matter  of  fact, 
the  second  variation,  A gCp,  of  the  drag  coefficient,  Cp,  resulting  from  an  ar- 
bitrary variation  in  7f  , of  Aft  , may  be  written  in  the  form 

/t  rf  rt 

/ w(u-$oa  ? ?C+)  Arf(^)d<a  cfy>0 

Jo  Jo 

by  reference  onoe  again  to  Eq.  (1).  (W) 

Now  if  the  conditions  stated  as  Eqs.  (6)  and  (7)  are  still  to  be  satisfied, 
then  A gCp  corresponds  to  the  drag  ooeffioient  that  is  produoed  by  an  annular 
duot  having  a value  for  R0  that  is  the  same  as  the  inlet  seotion  and  the  same  aa 
the  outlet  seotion,  aa  well,  and  whose  meridional  line  is  desoribed  by  moans  of 
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the  funotion  A >7  (z).  The  drag  coefficient  of  this  duct,  however,  is  always 
going  to  be  positive. 


4.  Solution  of  the  Variational  Equation,  Eg.  (8) 

The  solution  of  Eq.  (8)  is  carried  out  by  replacing  the  indioated  funotion 
F(z)  whioh  appears  in  this  Eq.  (8)  by  its  longer  expression  in  terns  of  Light- 
hill's  W-function  obtained  from  Eq.  (5).  When  this  is  done,  one  has  that 

where  the  taro  new  constants  standing  on  the  right  hand  side  of  this  equation 
are  still  being  represented  by  means  of  the  symbols  X ^ and  X 2 that  had  for- 
merly been  employed  in  Eq.  (8). 

The  expression  now  obtained  as  Eq.  (11)  constitutes  a Fredholm  integral 
equation.  Its  solution  may  be  written  in  the  standard  form 


W C /)  ri  C*A)cf<A  + A,  + X^  ^ 


(U) 


V - x,  ^ X £ 


c ^ ^ ^ 2 ^ i ct 


(12) 


where  the  resolvent  kernel  P (a, a)  is  defined  in  terms  of  a Fredholm  series, 
whioh  is  uniformly  convergent  regardless  of  the  value  of  X.  . It  will  normally 
be  true,  however,  that,  in  those  oases  of  praotioal  interest  in  the  solution  of 
aotual  problems  suoh  as  are  under  examination  here,  the  values  of  & will  be 
suoh  that  one  may  replaoe  the  Fredholm  series  by  a Neumann  series  and  reap  the 
no  mean  benefit  of  sharply  reduoing  the  tediousness  of  the  involved  oaloulatione 
whioh  would  otherwise  result. 

Thus  it  will  be  to  advantage  to  express  the  resolvent  kernel  P (z,s)  in 
terms  of  the  series 


(Z) 

k + 


£' 


(n) 

X (13) 
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wh*”  K.  (y,*)  - MO*-  yO 

K ^ C^jj)  — f K 1 1)  K ( t ) ) d t 

ra  L U*) 

k (%,*)  = F k r?,t)  K<h-°  (t^)dt 

In  order  for  the  Neumann  series,  Eq.  (13),  to  be  uniformly  convergent,  it 
is  sufficient  to  stipulate,  on  the  basis  of  the  condition  given  by  E.  Sohmidt 
in  Reference  3,  that  the  relation 


1< 


1 


/ / /Ywru~  yoj  lud. 


(15) 


should  hold. 

Now  in  the  present  case  it  is  true  that 


AV  ( | ^ \ ) < 0.5  + — -J~  — — I ~~  y I (16) 

where  W(  £ ) has  been  denoted  by  A.  If  one  remembers  that  for  £ > 4 the  corres- 
ponding value  of  A will  be  so  small  that  it  oan  be  neglected,  then  this  value, 
appearing  on  the  right  hand  side  of  Bq.  16,  which  is  larger  than  W(|s-s|  ),  may 
be  substituted  in  plaoe  of  W(|s-z|  ) in  Bq.  (15),  so  that,  in  fact  the  value  of 
the  double  integral  appearing  in  Bq.  (15)  will  be  bounded  by  the  quantity £ ^/Q  ; 


that  is. 


W (U-y\)J  d*dy  < < -j 


and  thus  it  follows  that  the  condition  (15)  may  be  reoast  into  the  simple  state- 
ment that  the  Neumann  series,  Bq.  (13),  will  be  guaranteed  to  oonverge  provided 


that 


£<  2 VT  - 5.656. 


(17) 
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Of  course,  the  radius  of  convergence  of  the  Neumann  series,  Bq.  (13),  is 
somewhat  larger  in  aotuality  than  is  allowed  on  the  basis  of  the  limiting  value 
given  in  Eq.  (17).  For  large  values  of  JL  , however,  it  will  be  more  convenient 
to  look  for  a solution  of  Bq.  (11)  in  finite  terms;  this  oan  be  done  by  replaoing 
the  kernel  W(  |s-z | ) by  a trigonometric  series  having  a finite  number  of  terms. 

If  the  series  is  to  have  m*l  terms,  then  at  each  of  m+1  seleoted  points  along 
the  | s-i  |-  axis  lying  within  the  interval  from  0 to  4 this  series  will  be  re- 
quired to  take  on  the  same  value  as  the  W function  has  there.  Thus,  in  sym- 
bolic form,  the  kernel  will  be  expressed  as 
m 

'jD  m s 


w c U-i*  i;  = 


- 2^  *n  cos  r C 2 -.*>7 

7>-0  L 


(18) 


wherein  the  coefficients  are  determined  according  to  the  formula  of  Clairaut 

/ (4} 

De  la  Vallee  Poussin'  ' as 


7»-/ 


An  2L  w(^r)  cos 


sirt  n WJol+_fzJ_)Z  W_L*± 

777  + mn 


a • / 


and 
^ m 


for  n = 


TP)-/  / ' / ✓ 

2.  *0  / . 1 ''  rr  ^ m / m 

/t»/ 


M19) 


If  the  expression  for  W(|s-z|  ) given  by  means  of  Eq.  (18)  is  now  substi- 
tuted into  Eq.  (11),  and  if  the  auxiliary  definitions  are  made  that 


»*-V-  J *,-'T  > -t9  " -f„ 


(20) 


then  Eq.  (11)  is  transformed  into 


7{(9)»  X'-t  Z Xt£*&  + 4mJ'  y ' A^CcosnB  cos  n&,+3in  n6  sm  n 6j)  d&, 

rTt  (Hi) 

f ST"1  1 

£*9+4*1  7 An  Hr,  cos  n 6 +J  An  Kt j sin  n 0 


= \,  + ZX. 
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(21) 


whs re in  Hn  end  Kn  ere  oonstente  whioh  ere  defined  by  the  reletions 

Hn  * T U 7?  & i ) cos  n9,  e/9, 

Jo 

r 1* 

KnmJ  (9,)  sin  hG,  c/9 / 

In  order  to  evaluate  these  constants,  one  may  multiply  both  sides  of  the 
expression  for  (©)  , given  es  Bq.  (11*) » by  cos  n 9^6.6  ^ and  sin  n&  jd0  ± , 
respectively,  end  integrete  between  the  limits  of  0 end  rc  • For  convenience's 
sake,  when  oerrying  out  this  operation,  one  mey  melee  the  definition  thet 
Bn  - t*9)  cos  n 9 c/9 


end 


C„  ^ ( A , + 2\Z£*G)  sin  nO  c/6 

Jo 

end  thus  B0  - tc  f A / + A 2 A*  rt  ) 


(22) 


end 


Bn  ~ ~~z  ^ i ^ 


while  JC , 


= 0 

- — X + 2.  X / * — 

n A/  * Az * n 

» _ 2L  z X A* 
n A z * 


for 

n 

odd 

for 

n 

even 

for 

n 

odd 

for 

n 

even  . 

Once  these  understandings  heve  been  agreed  upon,  it  will  be  seen  thet  the 
H end  E oonstents  ere  related  es  follows t 

~ B,t  + £ * A i K'c( T+~si  ~~  st- 1 J 


B/z  H n 


t>i  Ki  * Ci  * for  r -0,2,4,“' 

end  i ■ 1,3,5, •• 


(23) 


bi  Hi  * B,  + A*  A*  ( <• -*-/»  - 6-/t  ) 
b*  K*  -C*  +**£?  Hi  (d* 


(23*) 


for  i « 1,3,5, ••• 
and  r - 2,4,6, ••• 
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wherein 


be  - 1 - rt  l*  A9  *\ 

6^  * 1 “ -J-  4*  A si  for  r / 0 (24) 

end  ^4  -l  - -j—  >£  * A i . I 

It  is  evident,  therefore , that  the  set  of  2nH-l  equations  in  m+1  unknowns  H 
and  E may  be  separated  into  two  independent  systems  of  equations,  one  of  whioh 
serves  to  determine  the  unknowns  Hr  and  E^  , while  the  other  allows  the  de- 
termination of  the  unknowns  and  Ej.. 

This  decomposition  into  groups  of  similarly  determined  quantities  may  be 
oarried  even  a step  farther,  inasmuch  as  the  seoond  of  the  Eqs*  (23)  may  be 
written  as 


where  another  simplification  has  been  made  by  setting 


J_ [_ 

i + st  si -l 


(26) 


and  thus  if  Eq.  (25)  is  substituted  into  the  first  of  the  two  relationships 
given  as  Eq*  (23)  it  follows  that 


- a*  * t*  Y .Si Ci  ^ Ai 4r 


L, 

for  t - 0,2,4,' 


(23") 


r ■ 0,2,4,«o« 
and  i - l,3,5,o**  . 

This  latter  relationship  now  permits  the  determination  of  all  the  H-values  for 
whioh  the  subscripts  are  even,  and  after  these  quantities  are  confuted,  then  the 
second  equation  in  the  set  given  as  Eq*  (23)  may  be  utilized  to  find  the  values 
of  all  the  E-values  with  odd  subscripts* 
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In  an  entirely  analogous  manner,  if  one  firat  makes  the  simplifying  defini- 
tion that 

b. 


i 


C/t  i + n, 


< - St 


(26*) 


it  will  be  found,  by  working  with  Eq.  (23'),  that  all  the  H- values  having  odd- 
numbered  subscripts  are  given  by  means  of  the  delation 

t;HC  - Bi*  jfT  ahC„  -ii±i-Smz  T Aj,  tuT  Afbl-  Ht 

*—i/t  bst  *—t/l  *— '(■  O/g 

for  t ■ l,3,5,o»« 

r ■ 0,2,4,«*» 
and  i M 1,3,5,*»*  , 

Once  these  quantities  have  been  determined,  then  all  the  K- values  having  even- 
numbered  subscripts  may  be  found  by  having  recourse  to  the  second  of  the  set  of 
relationships  given  as  Eq.  (23'). 

Since  A^  and  A 2 are  to  be  quite  arbitrary  constants,  then  it  is  worthy  of 

note  that  the  method  just  described  will  no  longer  be  applicable  in  that  case 

* * 

where  Jl  might  take  on  that  value  for  whioh  7t  Jt  ko  * 1 • If  this  were  to 

happen,  it  would  have  to  be  true  that  this  critical  value  of  X must  be  auoh  as 
to  satisfy  the  oondition 


kV  C^)  d^  = Z 


(27) 


where  it  is  taken  for  granted  that  W(z)  is  to  be  expressed  in  this  instance  by 
Eq.  (18)  and  where  it  then  follows  that  k^  is  to  be  evaluated  from  the  formula 

•Jo 

Now  it  merely  remains  to  note  that  there  is  no  value  of  A for  whioh  the 
oondition  stated  as  Eq.  (27)  will  be  satisfied;  this  faot  may  be  verified  by 
oheoking  through  the  values  of  Table  I,  whioh  is,  in  effect,  a definition  of  the 
funotion  W(s) . This  ncm-existenoe  of  a suitable  behavior  for  the  W(i)  funotion 
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in  this  speoial  oaie  is  all  the  more  reason,  therefore,  why  the  other  b^  and 
br  ooeffioients  will  not  drop  out  whatever  the  value  of  JL  may  be,  and  conse- 
quently one  may  conolude  that  the  method  being  proposed 'here  is,  in  truth,  always 
going  to  be  applicable. 

Finally,  it  is  to  be  remarked  that  Eqs.  (23)  and  (23*)  are  written  in  suoh 
a way  as  to  permit  the  determination  of  the  Hr,  set  of  constants  on  the  one 
hand  and  the  set  H^,  Kr  on  the  other  hand  by  means  of  an  iteration  process, 
provided  that  the  numerical  values  for  the  coefficients  multiplying  the  unknown 
constants  are  such  as  to  ensure  that  the  method  itself  will  converge. 


5.  Determination  of  the  Constants  A j and  A 2 

The  results  derived  in  the  preceding  Article  will  now  permit  one  to  write 
the  solution  for  (z)  in  the  form 


(28) 


where  and  Fg  are  known  functions. 

Thus,  one  is  led  to  the  conclusion  under  these  oiroumstanoes  that 

i f 2 r?  r 2 

* “ B • I Ftty  dp+Xz  / Fzty  c/fi  (29) 

J o Jo  Jo 

and,  therefore,  by  applying  the  condition  of  "closure*  given  as  Bq.  (3),  it 

follows  that 


/ 


f ^ ~ 0 

Jo  Jo 


(3*) 


and  meanwhile  it  is  also  true  that  the  perimetral  condition  originally  stated 
as  Eq.  (2)  will  result  in  the  further  restriction  that 

rt  r . re  r* 

Af  I d^  I I d op  I Ffc  d m C (2*) 

J a Jo  Jfi  JfS 
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These  conditions  (2*)  and  (3*)  constitute  the  definitive  requirements  whioh 
will  serve  to  determine  the  constants  A j and  A g uniquely.  After  their  eval- 
uation has  been  made  in  this  way,  then  there  will  be  no  further  trouble  in  ap- 
plying the  result  given  as  Eq.  (29)  to  derive  the  sought  description  of  the 
meridional  line  directly. 

6.  Case  of  an  Annular  Duot  with  Conical  Type  of  Central  Channel  (Ri  / RP) . 

Extension  of  Lighthill's  Method  to  Determination  of  the  Flow  About  a Conical 
Type  of  Duot,  Before  Optimisation 

First  of  all,  before  tackling  the  problem  of  determining  the  best  duct  shape 
in  this  case,  it  is  essential  to  give  prior  consideration  to  the  problem  of  how 
one  needs  to  go  about  finding  the  flow  created  over  the  outside  of  such  a con- 
ical sort  of  annular  duct  which  arises  from  action  of  a uniform  supersonic 
stream  that  is  assumed  to  be  impinging  upon  it,  under  the  simple  assumption  that 
Rl  / R2  j i.e.,  where  the  radius  of  the  circular  cross-section  at  the  duct  en- 
trance is  R^and  the  radius  of  the  ciroular  cross-section  at  the  duct  exit  is 
Rg,  and  where  these  two  cross-sections  are  substantially  different  in  size.  In 
order  to  carry  out  such  a calculation  a procedure  will  be  utilized  which  is 
quite  analogous  to  the  one  employed  by  Lighthill^)  in  his  study  of  the  similar 
problem  whioh  is  concerned  with  the  determination  of  the  flew  about  a nearly 
cylindrical  sort  of  annular  duct. 

To  this  end,  one  may  start  the  analysis  of  such  a flow  problem  by  studying 
the  disturbed  field  that  is  produced  by  a distribution  of  supersonio  sources 
looated  along  the  x-axis  with  local  strengths  per  unit  distance  of  U—  f(4)  • 

It  is  well-known  that  the  body's  geometric  properties  are  linked  to  the  souroe 
strengths  under  these  circumstances  by  application  of  the  boundary  conditions 
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whioh  are  the  sane  for  any  and  all  arbitrary  general  meridional  planes  through 


the  body  axis.  Thus  it  may  be  stated  that 


Cx-Sj)  f (i,)  d i, 

Rl 


where  f (Sj)  • • 


(30) 


Now  if  JS  is  used  to  denote  the  semi-vertex  angle  of  the  cone  which  forms 
the  basio  shape  upon  which  the  aotual  exterior  contour  is  to  be  built  up,  it 
will  be  quite  legitimate  in  the  present  situation  to  let 

R e * R i ( X - BR  f)  • 

If,  furthermore,  the  following  changes  in  notation  are  made  for  purposes 
of  easting  the  above  expression  into  more  tractable  forms 

6,  « BR.  j t 

and  f (4)  • h (t) 
then  the  integral  relationship  written  as  Eq.  (30)  becomes  converted  into  the  ex- 
pression 

L I jjjLxzll  till  eLi 

* t+A8?  + (32) 


or  it  is  even  simpler  when  written  as 


9 (%)•  — 2 

J [(i-y6B)d  -t][d+j3 3)}-t+lJ 


If  one  now  multiplies. both  sides  of  this  latter  equation  by, 


teg rates  between  the  limits  of  0 and  s 

h (t)cft 


V 


(*)  * 


, it  follows  that 


dt 

i s-s 


and  in- 


(<*  + J-t ) 


t__J  [n-^B)rt][cn-/3B)^z] 

i-AB 

(32«) 
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as  a result  of  interchanging  the  order  of  integration  with  respeot  to  t and  s« 
Upon  carrying  out  the  process  of  integration  with  respeot  to  z which  is 
indicated  on  the  right  hand  side  of  this  Eq.  (32')  it  is  found  that 

( l -sSB)* 


/ i * 

*■ M ~H 


W)i 


1 + 


AB 


1-J3B 


/ -Jdz  B z 


/ -£B 


ctt 


(33) 


z*n:  K'*>~ 

VZa **  E(f7+. 

t-jse  F f/x* 

/ +AB  yx+- 

(34) 


in  whioh  the  complete  elliptic  integrals  of  the  first  and  seoond  kinds  have  been 
denoted  by  the  symbols  F and  , respectively. 

On  the  basis  of  this  definition  for  K(x)  given  as  Eq.  (34)  it  is  seen 


that  K(0)-V  1 B . In  addition,  let  the  following  auxiliary  definitions  be 


made  for  convenience's  sake i 

~ 1-sBB 


h n-AB)]  * h*  U) 
h [tj  C/-AB)]  « A*  it,) 


(35) 


and 


L 


^ U)  - iT~Fe 


7T 


i 1-jQ6^ 

It  then  follows  simply,  by  taking  the  derivative  with  respect  to  s of  both 
sides  of  Eq.  (33),  that 

• ( t,  \ 

1 


UT  (A)  ~h  *U)  + / h *(t,)  yy- 


■s6B<a  y/y ZSTT, 

dK 


dt, 


(36) 


where  the  dot  has  the  usual  significance,  i.e.,  K(x)»  '<£jf ' . 

Now  it  so  happens  that  Eq.  (36)  is  a Volterra  type  of  integral  equation  of 
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the  second  kind,  and  its  solution  may  be  written  down  at  onoe  as 
fy*C+J  » t tr  (*)  - -j===g f p (*,  tf  ) w-  ( t f)  dtf 

Jo 


wherein  the  resolvent  kernel  is  given  by  the  following  series,  whioh  is  onoe 
again  uniformly  convergent  in  this  case; 


n (*,  t()  =h* 


fT+~3TT 


t))‘t TTTb  h*(3) (*> (3Q) 


wherein  the  several  H*  symbols  have  been  substituted  for  the  following  more 


complex  expressions^ 
H*  (■+)  tf)  — 


(i n^stf) 


7 )+yfi  B tf  i 1l+^Btt 

tj)  » H*  U>u.)  H*  (tx^o  du 


H*(3) 


and  so  on. 


Now,  further  consider  the  first  of  these  integral  relationships:  in  more 


extended  form,  one  has  that 


H*(z> 


(i md 


1+J6d+  -iJij6Bu 


and  if  a change  of  variable  is  made  in  whioh  one  lets 

u - t / 

1+s6Bt,  * U * 


where  the  symbol  y has  been  employed  momentarily  to  stand  for  the  relation 


In  an  entirely  analogous  way,  it  will  be  seen  that 


and  likewise  for  the  others. 


On  the  basis  of  these  evaluations,  it  follows,  therefore,  that  the  value  of 
the  resolvent  kernel,  P (s,ti)  , may  be  obtained  from  the  series 
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and  upon  Insertion  of  this  result  into  Eq.  (37)  it  is  found  finally  that 


7**/  -7 ”«>>*'  l57,) 


7.  Confutation  of  the  Pressures  and  Resultant  Forces  Acting  on  the 
Cone-Like  Body 

The  component  of  velocity,  lying  in  the  direction  of  the  x-axis,  -which  is 
due  to  the  disturbed  flow  about  the  body  in  question,  is  given  at  any  arbitrary 
general  point  on  the  surface  of  the  cone  (which  in  characterized  by  the  fact 
that  the  given  circular  entrance  and  exit  openings  of  the  annular  duct  under 
study  constitute  directrices  of  this  conical  surface)  by  means  of . the  usual  form- 

nl*  i di  _ f *~a*»  / (jzjje, 

}x  ~Ja  r (X  - 


/ 1ft)  dt 

s ( i + y-t)*-  - a+j3B %)* 


h*  <t.)  <ft, 


where  the  free-stream  velocity,  T7«e  » is  used  for  non-dimensionalizing  pur- 
poses. 

By  insertion  of  the  value  for  h*(tj),  obtainable  through  means  of  Eq. 
(37),  the  expression  for  this  x -component  of  the  perturbed  velocity  field  over 
the  oonioal  surfaoe  (as  oompared  with  the  free-stream  velocity)  is  obtained  as 


I 


where  the  symbol  T(s,  u)  employed  in  this  integral  stands  for  the  mixed  ex- 
pression 


Now  with  a slight  amount  of  further  manipulation  it  is  seen  that 


where,  as  usual,  the  dot  over  a function  serves  to  indioate  that  the  deriva- 
tive of  the  funotion  so  marked  has  been  taken  with  respeot  to  the  sole  para- 
meter upon  whioh  it  depends. 
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It  thus  results  th&t  the  sought  x-ooaponent  of  velocity  is  given  by  the 
relation 


i avion  i- 

/ VT  t / ZJjnJLl  v,y 

~ d X ~ rt  V f-ySB  l_  ^ ° J /”£  5a 

■] 


2 T (<$,  cl)  du 


c/u  r 

~ h 


y u - ^2 


(40) 


i !+  X3  B 

Now  the  second  term  on  the  right  hand  side  of  this  Eq.  (40)  may  be  further 
reduced  by  interchanging  the  order  of  integration,  so  that  in  its  rearranged 


form  it  becomes 

VT 


/ 


^ T_  (2j,  u)t/u- 


TC 


V /-><?  5 / dls§Tu 


u *+  , . 


Yu- 


d* 


Yz_  j 

n -i  I see 


'*■(*)</*  f9 


V a --<4-  VT+/TB 


LL 


■I 


*8  * 

g (<*)  Z d^ 


where  Z (*,s)  is  the  integral] 


z 


^ T.  C&JJ±1 

V l rn^&eh  tc  i 1 -yS 8 


-fz 


(41) 


On  the  basis  of  this  definition  it  is  obvious  that  one  has,  in  addition,  that 


z ~^‘77^3  P 


T C ^ j u) 

-m  r 1 s§  ~b  a 


d u 


and  thus  one  may  write  the  expression  for  the  sought  x -component  of  velocity  as 

(40*) 


- " 9*(°)  z +/  j*  (+)  * (yi-A) 

r 

z * [9*U)]' 
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Further  information  about  the  funotion  Z (***)  may  be  obtained  by  in- 
serting into  Bq.  (41)  the  expression  for  the  T(z,u)  funotion  given  as  Bq*  (39). 
When  this  is  done,  it  turns  out  that 


provided  that  one  momentarily  makes  use  of  the  new  variable,  defined  by  the 

rf  - u 

relation  T = r~  ~ to~  • 
f + yOD  u 


Now  also  let  the  following  new  definitions  be  set  up  for  simplicity*  s saket 


m = 


_Z 

7 C&B 


_(i 


•jSd  t+AB 


H&5*  /sb)( £ 


r / )/±t*d  ) 

( * * AbJ(  / -s6B  Z f-J3B~ 
~b)(i^S3  ” H-sdB  i‘/-ysz6x  "*) 


1 * -f&B  (i  _ T . 2 ^)/3'r  +-L  \ 

[/-S00  f+£B  /(  tt£W  £dj 

It  then  follows  that  the  expression  for  Z (z,a)  reduces  to 


i7&5\ 
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Then  if  one  notes  that 

YT  / _ 

it  i T~Wb 


. , , ¥ $ G 

2 * “ 


Z i TZ&JL 

l+yti  By 


n 

z 


t By 


(44) 


it  is  easily  perceived  that  the  x-component  of  the  velocity  in  the  disturbed 
field  may  be  rewritten,  from  the  formulation  obtained  as  Eq«  (40’),  into  the 
quite  compaot  statement 


Uco  #x  * 9*  (zh+ydBy  ~ [ 3 ** 

L c 

~ $ Cy)  ~J  5 <f*  C*)  d-** 

where  the  partial  derivative  of  2 with  respect  to  s is  denoted  by 


(45) 


£5  * 


(46) 


(47) 


Corresponding  to  this  result  for  determination  of  the  x-component  of  the  dis- 
turbed velooity  field,  the  following  expression  for  the  local  pressure  coeffi- 
cient holds : 

-ffihr  •9<9)-p (?,*■) 3*  *■+ 

= ^ (?>■*)  ? 

where  an  improvement  in  brevity  has  been  wrought  by  setting  Sf  (y } *+)& (l+jQB+)§  )» 

Consequently,  the  drag  coefficient  produced  by  this  sort  of  flow  aoting  on 
the  oonioal  duct  is  given  by 

x>*  rr*<  r*' 

Co-4-jrl  / 


- ^ r pySBp)  ? £ ; d9p 5 V- 


V)0+j8B+)cU  I 


(48) 
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where 


S*  - 3(z  ,s) 


while  5 • S(s,z) 


for  0 & b & z 


for  z < s < 


8.  Setting  Up  of  the  Equation  VThioh  Will  Determine  the  External  Shape  of 
the  Conioal-Like  Body  Giving  the  Minimum  Wave  Drag 

By  use  of  the  expression  (48)  just  derived,  it  will  be  possible  to  set  up 
a variational  equation  which  will  allow  the  determination  of  the  required  meri- 
dional line  7 (z)  . This  is  done  in  an  exaotly  analogous  way  to  what  was  done 
in  the  case  of  the  cylinder-like  duct  treated  in  Article  3.  The  sought  meri- 
dional line  is  the  one  which  will  produce  a cone-like  annular  duct  having  mini- 
mum drag,  and  which  will  be  subject  to  the  same  conditions  as  stated  in  Artiole  1 
in  the  form  of  the  perimetral  conditions  (a)  or  (b).  The  stipulation  stated 
as  restraining  condition  (a)  is  converted  into  mathematical  form  in  the  present 
instance  in  just  exactly  the  same  way  that  was  met  with  in  establishing  Eq.  (2), 


so  that  here  also  one  has  that 


'V  l 


7£  (^)  c/j,  = C 


must  be  satisfied,  while  the  formerly  stated  perimetral  condition  given  as 
Eq.  (2'),  having  to  do  with  the  volume  enclosed,  takes  on  an  altered  form  in 
this  present  case,  because  now  it  is  required  that 

r*t  r v 

I ( l+yQBnp)  f (+)  d*  » D - constant  (5C 


where  the  approximation 

Re  * R i s* 

Z 


“ +/9  Cx-6Rf)  - 


has  been  utilized  onoe  again  just  as  it  was  in  writing  Eq.  (32). 
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It  must  also  bs  true  that,  whatever  the  variation  in  the  oontour  of  the 


duot  happens  to  be,  the  further  restriotion  stated  as 

'6R,+  L 

dRe 

BR,  ** 


of  X **  R % — R j * R j ^ 8Aj 


must  be  invoked;  this  ^closure*  condition  on  the  contour  may  be  simplified  to 
read  simply  that 


f] 

V (?>  d? 


E • cons toh t 


(50* ) 


must  hold. 


When  an  arbitrary  variation  to  7£  , again  denoted  by  ^ , is  brought  about, 

a first  variation  in  the  drag  coefficient,  C^,  again  denoted  by  will  re- 

sult, and  it  will  have  the  form 


The  terms  appearing  within  the  square  brackets  in  this  expression  may  be 
combined  in  the  following  manner*. 

Firstly,  the  order  of  integration  may  be  interchanged  in  the  first  of  these 
two  integrals , to  give 


ft  rJf, 

( I + ftp ) / ( / + ^ 3*)  *5  * A ^ (•*)  o(x>>  * 


C / A fa)  cf*  j (f+jSBp)  {ft)  of-g 
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r*f  ♦ 

integral  / ( t+/6Bp)S  ^ ofrj 


where  one  must  remember  that,  in  the  evaluation  of  the 

it  will  be  neoessary  to  make  use  of  the  oonvention  previous ly  stated  as 
Eq.  (49);  i.e.,  it  is  agreed  that  S*  * S(s,z)  for  0&z  ^ s 

while  * S(z,s)  for  b<i  & 4 \ 


Upon  interchange  of  the  s and  i symbols,  it  then  follows  that 


r*>  pi, 

/ ( / +/SB-4)  A y U)  e/<*  / Cf-tyBB#)  5*  7 f^)  c/y  = 

Jo  j J0  j 

J C1-t/6Bj)A  ? (y)  c/gj  (1+j0d+)3*  7 U)  d* 


in  which  the  convention  must  now  be  observed  that 

S*  ■ S (z,s)  for  0£  s 6 z 

while  *3  (s,z)  for  z s JL 


Consequently,  it  then  turns  out  that 


provided  it  is  understood  that 

Fty 

On  the  other  hand,  on  the  basis  of  the  restraints  given  as  Eqs.  (2")  and  (50') 


// 


( t +^0*) 


S*  f 


(+) 


(52) 
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it  follows  that  the  variation  must  also  satisfy  the  conditions: 


■A 


t?  C *)  - O 


Jo 
and 


r 

J A?  - o 


Hi 


or  else  the  simplified  versions 

^ d 7?  C<p)  c/y  O and  j 4 7?  » O 

Jo  Jo 

must  be  invoked,  while,  if  the  perimetral  condition  on  the  volume  enclosed  be 


(53) 


tween  the  duct's  inner  and  outer  surface,  as  stated  by  Eq°  (50), is  to  be  the  one 

governing  the  variational  procedure,  it  will  be  required  that 

£, 


J 1 of*  - O 


* 0 ''O 

or*  by  interchanging  the  order  of  integration  once  again*  this  minimum  volume 
perimetral  condition  states  that 


f '(*+  Ayi ( ^ ** 


(53*) 


must  hold. 

If  the  first  variation  in  C^,  denoted  by  A C^,  is  now  to  be  made  to 
vanish  by  proper  selection  of  7?  (2)  , so  that  regardless  of  the  variation  in 
7?  , denoted  by  A 7?  (*)  £so  long  as  Eq,  (53)  is  satisfied^  it  must  be  true 
that 

iff#)  - +(l+/3D'g)  F (^)  + */  ♦ £ m O 

(54) 

is  to  hold,  while,  on  the  other  hand,  if  the  constraining  conditions  to  be 
satisfied  by  the  variation,  At?(z),  are  the  ones  stated  as  Eq,  (53')  together 
with  the  second  one  of  the  set  given  as  Eq,  (53),  it  must  be  true  that 
a 0 7?  (p)  - +(f  +ja Xfi-  '?)-  o (54*) 

holds , 
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9.  Solution  of  the  Integral  Equation  (54)  for  Eq.  (64* )3  Giving  the 
Contour  That  Will  Produoe  a Minimum  Drag 
If  the  expression  for  F(z)»  stated  as  Eq*  (52),  is  now  inserted  into 


(54*) 


Eq.  (54)  one  obtains  the  relationship  for  (z)  in  the  form 

- if  (•+**+)  A w *+ 

which  is  a Fredholm  integral  equation,  and  it  is  the  basic  determining  relation- 
ship which,  together  with  the  conditions  written  above  as  Eq.  (53) 

£or  Eq.  (53'  )J  , serves  to  define  the  sought  function  for  the  meridional  line 
producing  least  drag,  ^ (z) , and  which  also  determines  the  constants  X ^ and 

X 2 (or  A 2 )• 

The  solution  of  this  Fredholm  integral  equation  may  once  more  be  represented. 


just  as  was  done  in  Article  3,  by  means  of  the  following  expression 

X / + X 


/J/ 

p(9*4  rhra 


ate 


(55) 


wherein  the  resolvent  kernel,  denoted  by  P (z,s),  is  to  be  handled  in  just  the 
same  way  that  the  analogous  situation  was  attacked  when  the  solution  to  Eq.  11, 
given  as  Eq.  (12),  was  encountered  in  identical  form  in  Artiole  4. 

A solution  to  Eq.  (54*)  may  actually  also  be  obtained  directly  by  pro- 
ceeding aocording  to  the  method  propounded  by  Goursat^®).  This  procedure  will 
afford  a treatment  of  the  present  problem  which  will  be  reoognized  as  being 
quite  similar  to  the  one  employed  previously  in  Artiole  4. 

First  of  all,  let 

(!  + /3B^)^  te)  - f(y)  and  ( <+)  =?(*) 

and  then  rewrite  Eq.  (54*)  in  the  form 


1 


*1 


• 5 * + x,  + % 


(56) 


How,  in  addition,  sot  0 - Z.  I and  9\  • % l and  let  S**be  represented 

* 1 ‘'l 

(approximately)  by  means  of  the  trigonometric  series  of  a finite  number 


(2m  + 1)  of  terms  in  the  form 

to 


?**  - 


to  ” 

7 A n cos  n 6 cos  n&i  + J A m + 

4t£  / 


^ Sin  k 9 Sin  kOf 


(57) 


where  the  constants  and  are  determined  by  applying  the  condition 

that  the  sum  of  the  squares  of  the  deviations  between  S and  S , through- 
out the  interval  from  0 to  7 1 , is  to  be  a minimum;  i.e<>,  the  quantity  E* 

is  to  be  a minimum,  where 

r 


(57«) 


is  defined  by 


Jo  J O 

and  where  the  difference  between  S and  5 

K,(#,9t)  - 3 *-  S**> 

If  the  expression  for  S**  given  as  Eq.  (57)  is  inserted  into  the  above 
expression  for  K » the  resulting  expression  for  E?  will  be  seen  to 

be  a function  of  the  A' s only,  and  the  condition  that  E*  should  be  a minimum 


(or  maximum)  is  given  by  the  statement  that 

<95* 


<9  A 


= O 


(57*) 


must  hold  true  4 
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Mow  let  the  following  convenient  notation  be  introduced: 


■ Tf 


rr  r*  K 

S* cos  n O cos  ft  Ot  d 9 d9  f 


for 


n « 0,1,2,* 


nit 

- / / 3* sin  k& sift  kdj  dO  d&f 

Jo  Jo 


for  k ■ 1,2,3,* **,m 


pic  nrt 

= cos  ft 

da  do 


9 cos  n Of  cos  i 0 cos  i Of  dO  dO  f 


Tf 


= O for  n 4 if 

7 1 ftlt 

cos  n&  cos  ft  Of  sift  bO  sco k&f  d&  cfOp 


Jk  aV  T 


for  k*n  odd 

(and  thus  also  for 
k - n odd) 


for  k * n even 

(and  thus  also  for 
k - n even) 


Am  + k.,n 

Zif  k"J  =o  u k / J 
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Writing  out  the  conditions  for  an  extremal,  Eq«  (57*),  in  explioit  form, 
one  obtains  the  following  relationships,  where  m is  taken  to  be  even,  for 


instance : 

€ ° ” Oj  o A o 7/^  O)  m + 1 ^ m + 1 '+^'o3  m ^ wt3  + ’*  a m-f  ^ z nn  -/ 

€1  s ^1,1  htnfZ^fintz  + + Jr****  2«?  A 2 fin 

^ 3 =/66-tJ  2 ^2  m + f ^ m + 1 ’ "ty^LZjZfir>  - f ^ 2.  fin  - / 


1(59) 


™-',o  +s“m  + hZ  ^Z^t,'VUynrhyr>Ayfi^^)J  7»rf  ^finr, 

£rn+z  * >'66  n»-rz,i  rn  +.  2,3  <^5  * * * V*4?*  *2.,  w-;  ^ »■>  - f "t/^rntZj  mt  z ^Tfii-Z 


The  set  of  equations  just  written  as  Eq.  (59)  may  be  separated  into  two 
independent  systems  of  equations.  One  of  these  systems  contains  m+1  indiv- 
idual members  of  the  type  given  by 

m-i 

Jk  A 

6, 


?n  * s**  *9,  n ^ A * ^ S**  nt  vn  + ^ 


and 


“ i 


fin 


> (59*) 


€m-t-  k.  *V^  7»  + k)  vn+k-  ^ fin+k’t^^/u,rn*k)  n ^ fi  for  a m 0,2,°°°  ,m 

nm°  and  k * l,3,*°°,m-l 

The  other  independent  set  of  equations  contains  m individual  members  of 
the  type  given  by  ^ 

* i m S*"  L > L *"  X"  ijfin-bJ  A m+  j 

and 

+j  m*j  ^-fin+ j j,  i ^or  i * 1#3,  * * 4 ,m-l 

4 * f and  J - 2,4,*  ° ° ,m 

and,  what  is  more,  these  above  two  sets  of  equations  can  be  reduced  to  the  fol- 


(59") 
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lowing  two  independent  sets,  one  of  which  will  have  $ . 1 individual  members. 


and  is  given  by 


ho  » f 


^ h * / /* hj 


z 

k~l 


no- 


' no+k  


>n+k. 


ho 

i 


dtsslLJL-  ^ 

s^m+kjto+k. 


’] 

(60) 


for  n ■ 0,2,***,m 
* 0,2,**<>,m 
and  1c  ■ 1,3,  * * * ,m-l 

while  the  other  will  have  only individual  members,  and  is  given  by 

M 


m 


£<  iji  'O i ~h^^/u ij  mtj 
y-i 


r m-f 

f / f>  V 

/ yU7»+j)  7r»tJ  / 

L p.f 


l j£  I 
^ 1(60') 

for  i - 1,3,*** ,m-l 
i - 2,4,°**,m 


and  p * 1,3, <’**,iii-l  . 

It  is  worth  pointing  out  at  this  point  that  a very  significant  property 
that  is  characteristic  of  the  E's  may  be  seen  to  be  brought  out  by  the  array  of 
equations  given  above  as  Eqs.  (56);  that  is,  it  is  quite  easy  to  deduce  from 
them  that  the  computation  of  the  E's  may  be  performed  without  even  the  neces- 
sity of  having  to  know  what  S*  is  explicitly. 

As  a matter  of  fact,  in  order  to  give  some  elaboration  on  this  statement, 
let  attention  be  focussed  first  of  all  on  the  distribution  of  supersonic 
sources  which  are  to  be  located  along  the  x-axis  and  which  will  then  produce 

the  required  radial  component  of  velocity  given  by  the  formula 
TJ  / \ tt  co5  n 6 

($)  * u-  /T^sB  H 6 

at  eaoh  point,  located  in  the  interval  BR^  £ x £.  L*BR^  , and  lying  on  the  gen- 
eratrix of  the  cone  which  is  determined  by  the  entrance  and  exit  openings  of  the 
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annular  duct  under  consideration;.  Also,  let  $ represent  the  potential  whioh 
describes  this  flow.  Then  the  relationship  given  as  Bq»  (37')  may  be  used  to 
find  out  how  the  strength  of  the  source  distribution  must  vary  along  the  x-axis 
in  order  to  satisfy  the  imposed  boundary  condition  at  the  surface  of  the  duct, 
or,  on  the  other  hand,  one  may  employ  the  approximate  procedure  described  by 
von  Karman  and  Moore  in  Reference  7 to  effect  the  fame  end  more  simply* 


manner  as 


The  corresponding  drag  coefficient,  that  is  defined  in  the  conventional 

fCo).  * 


-zn ■ RUi 


/ 

J B/ 


_L 


ft  Rj- 


1 p tJ  * 

2 -'«•  uo° 


(61) 


may  consequently  be  determined  without  further  trouble „ 

On  the  other  hand,  on  the  basis  of  a comparison  of  Eq*  (58)  with  Eq»  (48), 
it  would  appear  that  this  drag  coefficient  must  also  be  given  by 


<c °)n  ~ it  [ 4«  L 


cos  n & 

% e 


dO  - Z 


tL  m 1 
tr*  *n  J ' 


(61* ) 


Of  course,  one  may  evaluate  the  first  term  within  the  square  brackets  by 
means  of  the  following  reduction  formulas 

fo  ir^Tire  de'ir« coi  4^/4 - f2nrr)  '«•  *) J 

in  4r  [A‘  (-r1  i~Zn  n)-^i  ( 7 


+ n * 
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where  a m^3B  ^p-,and  C i and di  are  the  sine-integral  funotion  and 
cosine-integral  funotion,  respectively.  Thus,  it  follows  that  Eq.  (61')  may 
be  used  to  obtain  the  values  of  € n directly. 

In  an  entirely  analogous  mtnner,  one  may  obtain  the  6,^  values.  Be- 
ginning with  the  statement  of  the  condition  on  the  radial  velocity  components 
produoed  by  the  distribution  of  supersonic  sources  plaoed  along  the  x-axis; 


i.e,,  using  the  relationship  that 

it  S_L_n_k& 

+ k ~ ~QQ  ±L  g 

Ttr 


Ue, 


one  may  obtain  the  value  of  (CD)m+k  directly  by  means  of  a formula  that  is 
analogous  to  Bq.  (61)  where  is  the  corresponding  drag  coefficient 

that  is  defined  in  the  usual  conventional  manner  in  Just  the  same  way  as  was  done 


previously.  This  drag  coefficient  is  thus  obtained  as 

a*  r 


( co)mt  its  * 


**  / 4jP.  I sut  k e - 


z /4 

rr  Jo  6 c/9~  "re*-  e » 


•*; 


(61") 


and,  in  this  case,  the  evaluation  of  the  first  term  may  be  performed  by  means 
of  the  extended  formula 

/ 


rn  si»z  ke 
do  Tt 


9 


/ 2 k. 

de  - ^ cos  — 


-k SI>  f [J>- Ji  (¥)] 


so  that  it  again  follows  that  Eq,  (61")  may  be  employed  to  give  the  sought 
values  of  6 . 

^ It  may  be  of  interes't  to  point  out  that  the  first  term  on  the  right  hand 
side  of  Bq.  (61*)  £ or  on  the  right  hand  side  of  Eq.  (61") J will  give  the  value 
(continued  on  next  page)  - 39  - 


Returning  now  to  the  main  trend  of  this  argument  it  will  be  assumed  that  the 
two  independent  systems  of  equations,  Bq.  (60)  and  Bq<>  ( 60’ ) , have  been  solved, 
and  that  the  values  of  and  have  been  calculated,  once  the  approximate 

£4  ifc 

value  S has  been  accepted  as  a suitable  representation  of  the  quantity  S • After 
this  has  been  done,  then  the  determining  equation  for  eliciting  the  required  meri- 
dional contour  may  be  converted  from  the  original  statement  given  above  as  Eq<> 

(56)  over  into  the  new  form  which  appears  now  ass 

f(e)-  + ~ *)f  >**(*,)  c/e, 

- k,+2Az  Jt*  6 *J>f  (I'Z/SB/*  9)^ jfj  An  H„  cos  n 6i^AMtk  Kta^ej 


where  the  new  symbols  appearing  herein  are  defined  as 

✓,*  = Jj_ 

2 it 

Hh  •*  r1*  ^(&t)  c.os  h 

p Tt 

Kk  p(9,)  Sih  ko,  ctd , 

with  obvious  analogy  to  the  definitions  made  earlier  in  Eq.  (21) < 


(63) 


of  (*D).  which  is  the  value  that  Cq  would  have  if,  in  the  calculation  of  this 
drag  coefficient,  one  were  to  employ  the  formula  which  gives  the  pressure  under 
the  condition  that  flow  is  assumed  to  be  two-dimensional.  In  other  words,  it 
turns  out  that  the  € -values  are  proportional  to  the  difference  between  the  value 
of  the  drag  coefficient  (Cp)  , obtained  in  the  manner  just  mentioned,  and  the 
drag  coefficient  whioh  would  result  from  a calculation  based  on  magnitudes  for 
the  pressures  which  would  be  obtained  by  working  with  the  correct  formula,  given 
as  Eq.  (47). 
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Now  let  the  following  useful  constants  Bn  and  Cn  be  defined  as  follows; 


it 


p 

B n * / (X,*  Z t*  X 2 e)  cos  nddO- -ZZf  \~~i  > for  n odd 


* z> 


- \,n+  \J^,nj£or  n - 0 


= O 


j for  n even 


and 


Cn  (x,+ zs*\e)  sc*  *ecfe~£  x,  + z *eJx£^}for 


n odd 


- z j?T  x r ~ j f or  n ®ven 

' Y) 


*■  o 


for  n * 0 


with  obvious  analogy  to  the  way  the  similar  quantities  were  defined  in  Eq.  (22). 

Consequently,  by  utilizing  these  definitions  it  is  seen  that  the  following 
relationships  between  the  H and  K values  are  obtained  on  the  basis  of  what 
was  previously  stated  as  Eq.  (62): 


■Gft.  Hy£  ~ B/f  + B/fl  A /Tt-f  l K(+  ZySB/,  1 J 

■Ci  K 'L  * Cl  + i f byy  i A nnt/t  Z^dff 

x:t  H{  « S"*  b ^ A,?, +4  K/t+Z/dBJ^  Lt 

-o’*.  <n  -C*  *-**Y.bUAn+i  Hi*  is66e* c, 

mrnrnm  [ w* 


where 


and 

while 


and 


-c*  tr 

- y6d  Jt,*Z  Tt 

-C/t  */  ~£" 

-y6Bei*Z  n 

*1  * / -^+  JI  ^ • 

-s&BJ?*Z  it 

■cA-'-Slf  -jSBrf-n 

^ IS)  1-St 

A - / -*T  ^ ^ v i rvAn  + i 

A ’ n ^ — 

i 

st  - L 

h - > 

/ 

b,x  ■ it*  l 

— ^5 

L«  »/  H 

. r , i , , 

1 

(i-s*)*' 

(66) 


*21  .A""j  «j[-  f -'^rJ  f 


A/ 

Zst 


ft  *- 


m+st 


'*[  ( i i-*-)z  + (L  -se)zj 

+ L;.  ^ ^ f-  ^/;  /V"  f 6jJ~  *< H/  £ 


Li  * * (/-*)*] 


+Lj  *~J  KJ  [-f'm0‘*J  “z  £>  ij ]-4„+t- Kt-  jf. 

~1L;  Kl[ t/t*)2  * (7^T)21 

^ ^ [-  * 0 Jt * f f* 


for  J ■ 0,1,2, "‘.a 
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in  whioh  the  prime  on  the  summation  symbol  is  for  the  purpose  of  pointing  out 
that  the  term  which  has  both  of  its  subscripts  the  same  is  to  be  dropped 
when  the  terms  are  summed. 

It  is  admittedly  true  that  the  calculation  of  the  H and  K values  by 
means  of  Eq.  (65)  will  be  somewhat  more  complex  than  in  the  corresponding 
case  encountered  when  dealing  with  Eqs.  (23)  and  (23* ),  but  nevertheless  it 
is  to  be  remarked  that  Eq.  (65)  is  admirably  suited  to  treatment  by  an  iter- 
ation process  for  the  purpose  of  finding  the  H and  K values,  and  this  iter- 
ative procedure  must  surely  converge  if  is  small  enough.  The  successive 
steps  in  such  a process  will  consist  of,  first,  finding  a preliminary  eval- 
uation for  H and  K £let  them  be  called  ^ and  KU>J  by  dropping  out 
the  L- terms  in  Eq.  (65) , and  then  proceeding  to  find  a second  approximation 
by  utilization  of  values  for  the  L1 s which  are  obtained  by  inserting  into 
Eq.  (66’)  the  approximate  values  and  in  place  of  the  actual  H 

and  K values.  Once  the  first  approximate  values  of  L are  obtained,  then 
Eq.  (65)  can  be  solved  once  again  to  give  the  next  more  acourate  evaluation 
for  the  H and  K values,  and  this  process  may  be  kept  up  until  the  re- 
quired accuracy  is  achieved. 


10.  Determination  of  the  Constants  \ l and  A 2 

In  order  to  obtain  the  values  of  the  constants  A ^ and  A 2 one  may  do 
just  what  was  done  previously  in  Article  5.  That  is  to  say,  the  solution  of 
Eq.  (54)  £or  of  Eq.  (54 i)J  is  written  in  the  form 

~ X/  F! 

wherein  and  F2  are  definitely  known  functions. 
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Then  here  again  one  may  deduce  from  the  above  statement  that 

/-j  * f * X/  r *■  ty  *9 

Jo  J0  Jo 


and  besides,  for  geometrical  consistency,  it  is  required  that 

/ Rz-R  "'**  n*' 

B R, 


— • Q — v — * 

o rfit  r st 

x' I F’ 


(68) 


while,  because  of  the  perimetral  condition  stated  as  Eq.  (2*),  it  must  be  true 
that 


F,U)  </++  X. 


- c . 


(69) 


A i 

j r / 1 *•/  «•  *•*  ' '"-l 
'O  U o 

Similarly,  it  follows,  because  of  the  perimetral  condition  stated  pre- 
viously as  Eq.  (50)  that  one  must  impose  the  condition  that 

A,  f 1 Ct  + /SB £ )<*#[%  fvrySB#)  f^z  W ~ D ( 70) 

do  Jo  Jo  Jo 

if,  this  time,  the  volume  enclosed  between  the  duct's  inner  and  outer  surfaces 
is  to  ba  made  a minimum.  Thus,  this  Eq.  (70),  together  with  Eq.  (68),  wherein 
the  constant  A g has  been  substituted  in  place  of  Ag  , will  afford  the  means 
of  determining  the  constants  under  consideration  in  this  alternative  situation. 
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PART  II 


1.  Direction  to  be  Taken  in  This  New  Approach 

The  problem  which  was  stated  and  investigated  in  Part  I can  also  be  solved 
by  making  use  of  a mode  of  attack  which  is  quite  similar  to  the  one  used  by  the 
author  when  he  previously  had  occasion  to  study  the  problem  of  determining  the 
solid  missile  shape  which  would  give  a minimum  drag,  the  solution  for  which  is 
given  in  Reference  2.  This  new  alternative  way  of  treating  the  problem  now 
under  consideration  differs  from  the  approach  used  in  Part  I essentially  be- 
cause of  the  fact  that  in  this  new  manner  of  handling  the  details  of  the  so- 
lution one  proceeds  to  find,  as  an  intermediate  step,  the  distribution  of 
supersonic  sources  which  will  produce  the  annular  body  contour  which  will  have 
minimum  drag,  rather  than  seeking  out  straight  off  what  this  shape  of  external 
contour  must  be. 

2.  Computation  of  the  External  Drag  Exerted  Upon  the  Cone-Like  Duct 

Examination  of  the  more  general  case  where  R^<Rg  is  to  be  undertaken 
right  from  the  start  in  this  instanoe,  because  the  solution  of  the  problem  of 
finding  the  shape  giving  least  drag  whioh  pertains  to  such  a cone-like  body 
will  not  present  any  more  essential  complications  than  will  be  found  to  arise 
in  connection  with  the  more  restricted  case  where  the  basic  shape  is  almost 
cylindrical. 

In  order  to  derive  the  formal  expression  for  the  external  drag  coeffi- 
cient, the  apparent  mass  conoept  will  be  employed  as  follows i Consider  the 

mass. of  fluid  enclosed  between  the  planes  AA  A'  and  BB  B'  vhich  are 

o o 

normal  to  the  x-axis  (see  Fig.  2),  and  which  oontain  the  entrance  and  exit 
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seotions  of  the  duot  whose  external  contour  Is  to  be  determined  so  as  to  give 
minimum  drag.  One  of  the  other  two  enclosing  surfaces  for  the  mass  under  con- 
sideration is  the  simple  boundary  marked  out  by  a oircular  cylindrical  surface, 
O ->  , which  is  to  have  an  arbitrary  radius,  denoted  by  R' , and  whose  axis  is  co- 
incident with  the  x-axis.  And,  finally,  the  last  enclosing  surface  is  the 
exterior  surface  of  the  rigid  duct  itself. 


Control  Surfaces  for  Applying  the  Momentum  Theory  to  Determination 
of  the  Drag  Coefficient  Corresponding  to  a Cone-Like  Annular  Duct 

Now  lei.  Qg  dt,  Q'dt,  and  dt  denote  the  component  in  the  direction 
of  the  x-axis  of  the  momentum  possessed  by  the  following  masses  of  the  fluid, 
respectively!  the  mass  flowing  out  of  the  section  BB*  per  infinitesimal  in- 
terval of  time,  dt;  the  mass  flowing  out  of  the  oylindrioal  surface,  , 
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In  this  sans  time  unit;  and  ths  mass  autaring  tha  saction  AA1  in  this  inter- 
val.  In  addition,  also  lat  Fp  ba  usad  to  reprasant  tha  drag  whioh  is  craatad 
in  tha  dist&nca  batwaan  tha  end-points  of  tha  duct,  A and  8,  and  lat  Fp  ba 
takan  to  raprasant  tha  component  along  tha  x-axis  of  tha  forca  whioh  is  dua  to 
the  pressures  acting  on  the  sections  AA*  and  BB’ . Then,  on  the  basis  of 
the  momentum  theory,  it  follows  that  the  above-defined  quantities  are  related 
by  the  formula 

Fo  ~fP  <3,. 


(i) 


If  one  makes  the  further  definitions  that  U and  U are  to  stand  for 

^ r 

the  axial  and  radial  components  of  the  local  velocity  vector,  "ft  , respectively, 
then  it  may  be  seen  that 


rxB  nR' 

<?V  Qz-Q,-ZitR 7 S>UXU*  </x  + Zrt  />RUf  JR-rc -R*X 
Jr*  Jr, 


(2) 


>x/i 

On  the  other  hand,  the  equation  of  continuity  requires  under  these  cir- 
cumstances that 


Y r*f-Rf,)+f*'°R  v,  a'R.r/z'  r 
* R a JxA 


B 


y°  c/x  - o . 


(3) 


Proceeding  further  with  the  reduction  of  Eq.  (2),  the  following  auxiliary  defin- 
itions will  also  be  found  helpful;  i.e«,  let 

C/>  - OU  i-  u 
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Then  it  follows  that 


s v-[  (“-B*  -^r^RcfR** 'aI u K»  </x  . (5) 

**  Fo  J *A 


Furthermore,  one  has  that 


r*# 

— Ff>  * / ^*^bo R c/  R 

Jr 1 


and  since 


-/>.  U*.  u-P„  JjL-/^  Hp  +SU  Mm 


R' 


RU%  c/R 


the  expression  for  may  be  put  in  the  form 

A ^ A ^ ' 

u*6**JR- 

and  by  combining  these  expressions  according  to  the  dictates  of  Eq°  (1),  it 
is  found  that  r\  X & 


/t  J 


r /^?R,R(Bz  u*  t U Z)  JR 

J*z 


Let  it  now  be  assumed  that  R'  * Rg.  That  this  is  a legitimate  simpli- 
fication will  be  readily  admitted  provided  the  thickness  (height  between  the 
annular  duct's  external  surface  and  the  basio  internal  conical  surface  to  which 
the  circular  sections  AAQ  and  BBQ  belong)  of  the  configuration  under  con- 
sideration is  small  enough.  Once  this  agreement  is  settled  upon,  it  follows 
immediately  that  Fq  reduces  to 


- yO  * rr*  f**  JL.  . 
fi~  **■  u-  / t/_ 

JkA 


cfx 
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and  thus  the  expression  for  the  drag  coefficient  turns  out  to  be  simply 

- Fn  2 Si  FX°<L  .f_K  \ 

° n *°JX  v- 

Z U XA 

Now  it  is  well  known  that  the  u and  Ur  components  of  the  flow  are 
given  by  the  formulae* 

u _ _f*-g**  £ ( t)sr&) 

/ i (*-e,)*-e*Ki 


and 


v* 


]_ 


K'3** 

■Ux-4)*-  - B*- 


where  the  notation  is  the  same  as  that  employed  in  the  development  given  as 
Part  I of  this  paper. 

Thus,  if  one  first  lets 

U - L (t-^B) 


for  convenience’s  sake,  then  it  follows  that 


r - * 

D n R% 


B R 2 p X * B p X ~ 8 

I iL,  I 


Be, 


In  order  to  put  this  expression  into  more  traotable  form,  let  the  inte- 
gration with  respect  to  x be  carried  out  first  of  all.  The  result  will  be 


it  R. 


'L‘  rB*z  c*  -4z)  </* 


JCx  -$  )±3*zi’(*-4J‘--  5**J 


'</ 


'Lz+BRjl 


f(€L)  I Yfx-?,)*-#*#*;  T>t'x-4z)‘-3*  /Z 

k,*BRz 
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Tt  ~£l 


f Llf  'fm  ^r,7-4,> 

Jo  lJ£,  Jf+BR, 


+ f &,)+€, 


L^BRz 


£ 

rv-  4±)cf*  „ i (6?) 

i fx -4  f)  a-  3^f‘  7?x - 4 J *-  0 4 * J / • 


The  very  last  integral  in  the  above  formulation  may  be  simplified  further 

by  proceeding  as  follows ; _ 

n L ofz^ 

Let  u \ ^ *. 


*Z  {*»**.) 


>£,  + B2z 


1 Tx  - ^f)z~8 'i (x -%)*•- B 


and  then 


Hz  (tf>tz) 


ft/* 


ft  *.-**)  *9*. 


1+  to 


(7) 


where  the  improvement  in  symmetry  has  been  made  by  applying  the  transform- 
ation defined  by 


Vz  8R2 
+ _ €*. 


and 


5/22 

''  -& 

K - — /,  r/->s£; . 


( 7 * ) 


The  evaluation  of  HgCt^tg)  may  be  performed  by  recourse  to  elliptic  inte- 
grals; i*e«,  one  finds  that 


* 
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Nw 


u ~T*J 


- tf)  i_7  * Vg  -r,  I j 


for  ^2-^1  ^ , 


1)  j kf)+  (tz-t f)  JJfPuk i)  7>f ) 


for  — 1 < i 
2 

where  IF  ■ F(JP0  ,k)  is  the  elliptic  integral  of  the  first  kind  with  argu- 
ment y>0  , which  is  defined  as 

y m wc,in  /S THU  Hi  ~ 15 

* Vr/jf - f,  * (tl  -t, i-z) 

and  the  modulus  of  this  elliptic  integral  is  denoted  by  k , where 

/ 2 \2 


— V ^ , 

The  other  symbol  employed  in  Eq.  (8),  i.e.,  y/  (^0,k,n)  » stands  for  the 
elliptio  integral  of  the  third  kind,  whose  corresponding  argument,  y ^ , and 
modulus,  k , are  given  by  the  same  formulae  just  set  down  above  as  Eqs.  (9) 
and  (9'),  while  the  parameter,  n , which  appears  in.  connection  with  this  in- 
tegral has  the  definition: 

"•-("rrr)  • 

***''  (9»; 

Finally,  the  other  elliptic  integral  of  the  third  kind  employed  in  Eq.  (8)  has 
an  argument,  modulus,  and  parameter  which  are  determined  aooording  to  the  re- 


lationships 


arcsin 


z(/l  - 
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In  an  analogous  manner,  the  other  Integral  appearing  In  Eq*  (6*)  may  be 
handled  by  letting 

^ / + no  - 

OlzJLlI  qLl 


Lz  + 


"t  ( tu  *t)  «, 


>4 1+3*2. 

and,  by  use  of  the  transformations  (7*),  it  then  results  that 


H i ( 1 1 > ?z)  ~~ 


tez  - 


/ f + t ■ , 


77^7T^rT='  7W*.  -*»)*■  -> 


and  the  evaluation  of  this  function  niay  also  be  effected  in  terms 

of  elliptic  integrals  as  follows: 

-2  ( 


»i  (*,.**)•  —t*  z *1  — Fft  ‘k)  + TT  »*) 

1 / *2  V *£  ' 


for 


ti-t2 


? 1 


- F(p3  , k,)+  ft,  -t2jW  (r3,k,,n3) 


(ii) 


for  tl“t2  < 1 
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If  one  now  scrutinizes  carefully  the  expressions  given  for  the  arguments 
and  parameters  in  Eqs.  (9),  (10),  and  (12),  it  will  be  seen  that  the  following 
interrelationships  exist  between  the  later  set  and  the  ones  given  previously! 
P Z.  > *&)  ■ Po  (*Z  J tf)  1 

y* 3 ( ^7  •>  ^z)  m ^ Z > ^ (13) 

»Z  (t/jtz)  -V  / 

(+tj+Z)  m »/  (+X*  t,)  . 


Consequently,  the  expression  for  the  drag  coefficient  given  as  Eq.  (6' ) reduces 


now  to 

Co  * 


(14) 


where  the  understanding  is  made  that  the  interval  of  integration  from  0 to 
^2  , both  in  the  case  where  the  integrand  is  f(tg)H(ti#tg)  as  well  as  in 

the  case  where  it  is  t (t^)H(t^,tg)  , will  have  to  be  divided  into  two  sub- 
intervals, which  will  run  from  0 to  t^  and  then  from  t^  to  in  the 

oase  of  the  first  expression  given  for  , while  in  the  seoond  expression, 
the  splitting  up  of  the  0 to  Q \ interval  will  be  done  at  the  point  tg  $ 
so  that  the  two  sub-intervals  run  in  this  oase  from  0 to  tg  then  from 

tg  to  Jl  2 • Of  course,  in  eaoh  oase  the  H function  to  be  used  will  depend 
on  whether  t^  is  less  or  greater  than  1 gj  i.e.,  H will  be  set  equal  to 
if  tg<  tj  , while  it  will  be  replaoed  by  the  Hg  expression  in  the  part 
of  the  interval  of  integration  where  tg>tj  • 
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3.  Determination  of  How  the  Supersonio  Souroe  Distribution  Must  Vary  Along 
the  x-jbcie  in  Order  to  Produce  a Minimum  Drag,  (C^)min 

The  statement  of  the  boundary  condition  for  the  external  flow  over  the  duct 
in  question  results  in  the  well-known  relationship 


/ 


*Es  « x 

cf  x R/+s<S  (x-6Rf) 


X- BR 


(*-$)  f (€)<*' 4 


Trt^7u~- 


where,  for  brevity's  sake,  the  upper  limit  on  the  integral  has  merely  been 
denoted  by  x - BR.  This  condition  may  be  improved  slightly  in  form  by  making 
the  substitution  CBJ  = Rj/kg  , to  give 


(15) 


and  consequently  a more  directly  useful  form  of  this  boundary  condition  re- 
sults as 

Rm-Ri  _-/  ^ eft.. 


BR 


±-  = 5 


ft,-a>X/-sS3)  . 

rt4-t*.)  f Cl  J 


z 


' ao 


CM>  +J3  8( tf-co)  I y 7t,  - ~)  ~ Lev  TfiW(  tra$]  * 

(16) 


If  the  order  of  integration  in  the  above  expression  is  inverted,  it  then  may 


be  written  as 


p r& % - *>)(> f & r #z 

RdR*1  mBJ 

The  next  stop  in  rendering  this  expression  into  more  suitable  form  for  later 
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manipulation  is  advanced  by  making  the  following  substitutions  t 

- u>) 


t*  ■ cu  + 

tr 

and 

*««• 

1 

ft 

where 

ii 

1 

— u>  • 


/SB 


- t,  m 


tL-  Y 


>86 


1-/36 

LO  • .i£r- 

/3B 


(17) 
(17* ) 


and  thus  the  integral  with  respect  to  t^  indicated  in  Eq.  (16')  may  be  recast 

into  the  form 

I 


f. 


Vz 


(t±zljA 


L 


[oo  tySB  (t,- oojjiftrt^)2-  [cjo-h/86(tl-ou')]z' 


1-^8 


•f  QO 


1_ 


ctt, 


7^86 


— r 


>3  BY 

/-/SB 


c/t 


u c t 

t*  1 !-y0xB*-J-Z  sSBYtf+yB*  B^Y*- 


* /,  * U 

..  /-y83 

where  ^ » ?z  - ce>  • —-£ 


(18) 


It  will  be  clear  that  the  evaluation  of  the  integrals  1^  and  Ig  re- 
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metral  condition  to  be  applied  below,  takes  the  form 

_ _ rfpt  ~ . 

- S [ IfSx, 


(19) 


■where  the  indicated  l(  z2»t2^  function  is  the  one  expressed  as  Eq.  (18' )• 
Now  the  desired  expression  for  the  integral  of  this  thickness  parameter 

R-Ri 

x is  given  in  extended  form  as 
BR2 

pBR,  + L pcv+fz  n z ~ &>)(* ~AB) 

/ ^ <*-«***/  dVi  / 1 (?#**)*<**')  **>■ 

d SR,  * d co  Uo 

p£* 

= ^Ctjf  (tt)dtz 


where 


X(tz) 


CO  + 


t-^SB 


1 tz)  2.  . 


I-JBB 


+ CJO 


(20) 


1 


Finally,  therefore,  the  statement  of  the  (a)-type  of  perimetral  condition 
discussed  in  Article  1 of  Part  I may  be  put  into  the  single  form 

f* 

■&  (t  z)  f (tz)  eft  z.~C  • (21) 

•O 

In  an  entirely  analogous  manner  the  other  "volume*  perimetral  condition 

may  be  given  a simple  formulation  by  first  noting  that 

pBRji-L  rBRj+L 

* J BRf 

where  onoe  again  the  average  duct  radius  is  expressed  as 

*«**<’  - R,  + /3(x-J*,) 


«■  56  “ 


and  the  internal  radius  is  given  in  terms  of  the  mouth  radius  by  the  relation 
* R-f  + A & • 

Consequently,  the  (b)-type  of  perimetral  condition  discussed  in  Artiole  1 
of  Part  I may  be  expressed  in  the  form 

r^'-kr.r  , 1 . r'?*~eu¥{7e,f\ 


p ootQO  r 1 P & ‘ 

/ B**la »(1-J38)+£B9A  lC#z,t2)f 

°°  Jf*  “0 

« r *-u* ' tz)f  ( tt ) dt * - d 


(*t)  dtt 


provided  it  is  understood  that  G*(t2)  stands  for 

. ££.  _ + ca> 

n 1~y<3B  r 


• j f Jl(#2  > tz) 

J £jt .*  co 


. (22*) 


t-jOB 

Of  course,  the  final  ^closure*  condition  to  be  satisfied  in  all  events  is 


Rj i-  & 

BRZ 


L-  - j-  (I-  *>)‘B  f . <*>-.  tz)i(tz)  d*t  ’ 00"t“t- 


The  first  variation  in  the  drag  coefficient  which  results  from  a var- 

iation  in  f (t2)  . denoted  by  A f(t2).  will  be  seen,  by  reference  to  Iq.  (14), 
to  be  expressed  by  the  following  equation 

ACd  - -ffc/ j'fA-r(t2)dtz J 1fCt2jH(t,3ti)  dtt 

* J , f J dt, J 


Z B±_  Rjl 

re  R.% 


- 67 


By  remembering  how  the  expressions  for  and  H g were  given  further  elab- 

oration  in  Eqs.  (7)#  (8),  and  so  on,  it  is  of  aid  for  the  present  purposes  to 


realize  that 

In  the  case  where  tj^tg  S : 


H ft  i jt^)  H ( } tf)  — H % 

_ tiAihzJjl 

- *1 


(tfj  1 1 i)  + H f (f  2 ) t/ ) 

F (?„,  k)t  : 4*  TTf-P^k.n) 

TZ  zt 


for 


t2-tx 


> 1 


- *)F(r,Jk,)i-z/'ti-tf)'nf?>,*l,rl) 


for 


t2-ti  ? i 


In  the  oase  where  0 < t2<  ti  : 


>(26) 


*%-£*  V (*•*>*>) 


for  tl't2 


> 1 


(t,  - i-2  - « )FfJ%^,)  +z(t, 


for 


Recognizing  that  the  H(t^,t2)+H(t2»t^)  function  just  discussed  will  take 

on  the  alternative  expressions  given  above  in  the  separate  parts  of  the  interval 

from  0 to  A 2 > let  the  convention  now  be  i&ade  that  this  function  is  to  be 

# 

denoted  more  simply  by  the  single  symbol  V(t^,t2) 
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Consequently,  the  first  variation  in  the  drag  ooeffioient  may  be  written 


as 


acd 


ZB1 

re 


o«  a* 

**■  r*z 


o 


(27) 


provided  that  the  £(tg)  function  really  defines  the  way  the  distribution  of 
supersonic  sources  must  vary  along  the  x-axis  in  order  for  the  corresponding 
drag  to  be  a minimum,  inhere  this  minimum  is  denoted  by  the  symbol 

Sinoe  it  is  required  that  the  conditions  expressed  as  Eq.  (21)^or  else  by 
Eq.  (22 )J  in  conjunction  with  Eq.  (24)  are  to  be  satisfied  at  all  times,  it 
follows  that  it  will  have  to  be  true  that 

Af  (t,)  I (£**  O)  dt,  - 0 


and 


Af(tf)  Mftf)  c/tf  » O 


must  hold  for  any  suoh  variation  ^ f(t^). 

It  will  have  to  be  true,  consequently,  that  the  function  f(tg)  must 


satisfy  the  relationship 

9*  * 

(rn  r2) 


I 


(28) 


On  the  other  hand,  if  it  is  condition  (22),  instead  of  condition  (21),  which  is 
required  to  be  satisfied,  an  equation  exactly  similar  to  this  latter  Eq.  (28) 
will  have  to  hold,  where  the  new  condition  on  £(t£)  will  be  obtained  from  the 
previous  one  by  the  mere  substitution  of  the  function  <J*(ti)  in  place  of  the 
function  fl(tjj  that  appears  here  in  Eq.  (28). 
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4.  Solution  of  the  Integral  Equation,  Eg.  (28) » 

The  determining  equation  for  the  required  eouroe  distribution  whioh  has  now 
been  derived  as  Eq«  (28)  is  an  integral  equation  of  the  first  kind,  and,  as  is 
well  known,  suoh  an  equation  does  not  have  a solution  in  general,  if  the  given 
function  is  arbitrary.  The  investigation  carried  out  in  Part  I of  this  paper, 
however,  showed  that  there  does  actually  exist  a function  which  describes  the 
external  contour  of  the  duct  in  such  wise  that  the  drag  produced  by  this  duct 
shape  is  a minimum.  In  addition  to  acknowledgment  of  this  fact,  it  may  also  be 
observed  that  to  each  such  external  contour,  described  by  any  such  meridional 
line,  there  corresponds  one  and  but  one  single  distribution  of  supersonic  source 
strengths  along  the  x-axis,  denoted  by  f (t)  , as  may  be  seen  by  referring  back 
to  what  has  already  been  deduced  in  Article  6 of  Part  I,  Upon  realization  of 
the  import  of  these  facts,  therefore,  it  may  be  concluded,  even  without  the 
necessity  of  carrying  out  the  proof  starting  from  first  principles,  which  would 
be  rather  laborious,  that  there  does  exist  a solution  to  Eq.  (28)  and  that  it 
will  be  unique. 

Before  going  on  to  the  determination  of  the  solution,  which  one  is  assured 
exists  uniquely  in  light  of  what  has  just  been  said,  it  is  worth  pausing  to  ob- 
serve that  the  kernel  of  the  integral  equation,  Mft^tg)  , exhibits  a loga- 
rithmic singularity  at  the  point  where  t^  • . In  more  precise  analytic  terms  it 

t 

may  be  stated  that 

- to?\s\  o(D  + <y(s)  (29) 

where  the  symbol  |cf  | is  used  to  denote  the  absolute  difference  between  the 
t-values,  i*e.,  i ill 

K-*>- k • 
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Let  the  integral  equation  under  consideration  be  rewritten,  therefore,  in 


the  form 


/ *(*>■)  \t[  ^|~  iffZgg  + *(*,) 

^ (28*) 


Is  a result  of  this  manipulation,  therefore,  the  funotion 


M(ti,t2) 


remains  finite  throughout  the  whole  interval  running  from  0 to  i jj  » and  it 

takes  on  the  value  unity  at  the  place  where  t£  ■ tj. 

In  order  to  obtain  an  approximate  solution  to  Eq.  (28')  one  may  prooeed 

by  first  dividing  up  the  interval  from  0 to  2 g into  m sub- intervals  of 

equal  length  to  be  represented  by  the  symbol  a.  In  addition,  the  function 

**!;*»>  is  to  be  replaoed  by  a step-function  (holding  for  any  arbitrary 
lag  |t2-ti| 

general  value  of  tj  whatsoever)  which  may  be  set  up  as  follows; 

U) 

Let  t^  represent  the  value  of  t-i  which  corresponds  to  the  midpoint 

(j) 

of  the  i-th  sub- interval,  and  let  tg  represent  the  value  of  tg  vhich 
corresponds  to  the  midpoint  of  the  j-th  sub-interval.  Then  the  definition 


is  made  that 


lof  t,  -riJJ 


Now  in  each  one  of  the  m small  sub-intervals  the  value  of  f (t  ) ia 

2 

taken  to  be  oonstant  and  to  have  the  value  which  this  £(tg)  funotion  attains 
at  the  midpoint  of  the  smajLl  sub-interval  under  examination.  If  one  then  rep- 


resents these  oonstants  in  the  form 


aj 


it  follows  that  the  integral  appearing  in  Eq.  (28*)  may  be  represented  by  the 
following  sumt 
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f)  ru)  t.9 


Ci)  | , 

*/  ~*z  I dtz 


J=/ 

where  ft  (ifO}  t£/)j 


'fi  AJN(t“>>  **J>J  *(+tllrJjy)+i  a[‘°?(i. 


t .2 

l Z 


. (j)  Q_ 
z ~ z 


to? 


f _ *■  fO 

Tz  ~z 


-n 


tP-t'0 


\t“Lt'0  |*f) 


(32) 

for  j / i 0 

If  one  now  inposes  the  oondition  that  Eq«  (28')  must  be  satisfied  at  all 
m midpoints  of  the  individual  m sub-intervals  just  defined  above,  the  fol- 
lowing system  of  m equations  in  the  m unknowns,  A , will  thus  results 

T N (tft  tij))  A (t'°,  t[i})  A]  t [uf  f - ,]  A; 

for  i • 1,2, ••*,»• 

In  con  junction  with  the  condition  (33)  it  is  also  necessary  to  have  the 


(3 


following  two  conditions  satisfied: 


constant  • 


T l( yr£d  + *(ZJ , ) a 4J  * constant  I 

fe  ^ r (34) 

and  **(*1^)  **Aj  m constant. ,J 

«/»/ 

Thus  Eqs.  (33)  and  (34)  constitute  a system  of  m+2  equations  in  the 
m+2  unknowns  A,  X x » and  X 2*  811(1  thus  th®y  afford  the  1988118  of  8i*rl-vin6 
at  the  solution  to  the  problem  in  hand. 

It  should  be  pointed  out  at  this  juncture  that  on  the  basis  of  the  defin- 
ition of  the  function  M(tlPt2)  that  is  given  as  Eq.  (26)  above,  and  if  one 
also  recalls  how  the  values  for  jp Q and  y>  2 were  originally  set  up  by  means 
of  Eqs.  (9)  and  (12),  it  will  be  readily  appreciated  that  when  the  i-subsoript 
and  j-subscript  are  widely  different,  then  the  value  of  M(tj  ,t2  ) will 

| (i)  ( j)  | 

lie  close  to  zero  for  large  values  of  | tj  - t2  | • Consequently,  it  then 

becomes  abundantly  clear  that  some  of  the  coefficients  contained  in  the  set 
A,  appearing  in  Eq.  (33)  overshadow  those  for  which  the  i-subscripts  are  close 
to  the  j-subscript.  This  interesting  property,  which  is  an  essential  charac- 
teristic of  the  Aj- values,  will  stand  in  good  stead  when  it  comes  to  solving 
the  system  of  equations  written  out  explicitly  above  as  Eq.  (33),  because  the 
method  of  iteration  which  would  be  employed  in  carrying  out  this  sort  of  com- 
putation may  thus  be  short-cirouited  or  telescoped  in  great  degree. 

One  should  not  fail  to  observe  the  very  significant  faot,  finally,  that 
the  parameter  n,  whioh  is  one  of  the  factors  whioh  identifies  and  specifies 
oonoretely  the  elliptio  integral  of  the  third  kind  that  enters  into  the  def- 
inition of  the  M function  under  consideration,  always  takes  on  values  that 
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are  less  than  -1  in  the  present  instanoe.  The  importance  of  this  statement 
lies  in  the  realisation  that  one  may  thus  proceed  to  evaluate  this  elliptio  inte- 
gral in  an  especially  easy  way,  i.e.,  by  utilization  of  the  simple  equivalent 
formula  established  by  Huel^.  In  order  to  evaluate  this  elliptic  integral 
under  the  present  circumstances,  therefore,  it  is  merely  necessary  to  obtain, 
as  a preliminary  step,  the  values  of  OL  and  m from  the  following  formula- 
tions s f 

71  -<Anz  (&  j h) 

and  = cj2_C£LxMll  gCP-IsZL  ihj  (35) 

__  - -T-  -gj 

where  sn,  cn  , and  dn  are  the  Jacobi  elliptic  functions  known  more  precisely 
as  the  sine-amplitude,  cosine-amplitude,  and  delta-amplitude  functions.  Fur- 
thermore, it  is  then  useful  to  define  quantities  q and  AL  according  to 
the  scheme 

/ OL 

9 \ F 

and  _j\  m.  ^ i kl 

zf 

so  that  finally  the  evaluation  of  UT  (JP, n,k)  is  obtained  from  the  relation* 


m Tf  (%  ri  j kj  « ^ toy 


Vj.  (AkLx) 

iTj  (A-  9) 


*4  (9) 


A 


(35') 


where  I?  denotes  the  Jacobi  theta  funotion,  and,  of  course,  as  previously,  F 
represents  the  complete  elliptic  integral  of  the  first  kind  with  modulus  k , 
while  F(>\lc)  is  the  elliptio  integral  of  the  first  kind  with  argument  A 
and  modulus  k . 

Fortunately,  therefore,  the  evaluation  for  the  elliptio  integral  of  the 
third  kind  appearing  in  Eq*  (26)  may  thus  be  carried  out  by  means  of  Eq.  (35') 
merely  through  reference  to  standard  tables. 
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5.  Determination  of  the  Meridional  Profile  for  the  Annular  Duct  Tfhioh 


Will  Have  the  Minimum  Amount  of  Dr* 


Once  the  distribution  of  supersonic  sources  located  along  the  x-axis  has 
been  obtained  as  explained  above,  then  the  determination  of  the  shape  which  the 
corresponding  external  duct  oontour  will  assume,  to  produce  least  drag  under 
the  stipulated  conditions,  will  require  but  a slight  amount  of  further  treat- 
ment of  the  information  already  at  hand.  As  a matter  of  fact,  all  that  is 


necessary  is  to  make  use  of  the  expression  given  previously  as  Eq*  (19)  for 

R -R 

the  thickness  parameter,  6 1 

br2 


v 

* 
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